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Strong coupling between electronic and mechanical degrees of freedom is a basic requirement for 
the operation of any nanoelectromechanical device. In this Review we consider such devices and in 
particular investigate the properties of small tunnel-junction nanostructures that contain a movable 
element in the form of a suspended nanowire. In these systems, electrical current and charge can be 
concentrated to small spatial volumes resulting in strong coupling between the mechanics and the 
charge transport. As a result, a variety of mesoscopic phenomena appear, which can be used for 
the transduction of electrical currents into mechanical operation. Here we will in particular consider 
nanoelectromechanical dynamics far from equilibrium and the effect of quantum coherence in both 
the electronic and mechanical degrees of freedom in the context of both normal and superconducting 
nanostructures. 
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I. INTRODUCTION 

One of the main goals of contemporary nanophysics is to 
attain the means for controlled operation of mechanical 
devices on the nanometer length scale. The required cou- 
pling between electronic and mechanical degrees of free- 
dom can be dramatically enlarged in conducting nano- 
electromechanical systems (NEMS) where electrical cur- 
rents and accumulated charge can be concentrated to 
small spatial volumes. Of particular interest in this Re- 
view are small tunnel structures that incorporate a mov- 
able element in the form of a suspended nanowire. This is 
because here, as we shall demonstrate below, one is able 
to harvest a host of mesoscopic phenomena for the pur- 
pose of transduction of electrical currents into nanome- 
chanical operations. 

One example of a nanoelectromechanical (NEM) de- 
vice is the nanoelectromechanical single-electron tunnel- 
ing (NEM-SET) transistor. Such a device may be built 
by self-assembly, employing mechanically soft organic 
molecules anchored at metal electrodes while attached 
to a small conducting "Coulomb dot" or "single-electron 
box'f^. In this way mechanical displacements of the dot 
with respect to the conducting leads are possible (see 
also Ref. [3j). Other approaches to building NEM-SET 
devices involve the use of a flexible nanopillar placed be- 
tween a source and drain electrode^ and placing a metal 
"dot" on the tip of an externally driven cantilevei^l Gen- 
erally, in NEM-SET devices the transport of electrons is 
due to a combination of single-electron tunneling events 
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(between either the source or the drain lead and the dot) 
and the mechanical transportation of charge through to 
the motion of the dot. 

Mechanical vibrations can be coupled to the electron 
transport both via uncompensated charge (electrostatic 
coupling) and currents (magnetomotive coupling). An 
electrostatic coupling can be realized in devices where 
a nanomechanical resonator is used as a movable gate 
electrod d^ l ^ l ^ i In such nanoelectromechanical devices, 
where the movable part is effectively "zero-dimensional", 
one has found that the electromechanical coupling can 
— under certain conditions — have much more dramatic 
results than simply vibron assisted tunneling, which fol- 
lows from a straightforward application of perturbation 
theory. In particular, it has been shown that an electro- 
static coupling of mechanical and electronic degrees of 
freedom can induce an electromechanical instability, re- 
sulting in self-sustained temporal oscillations in both the 
mechanical and electronic characteristics of the system. 
This instability is behind the new "shuttle" mechanism for 
mechanically promoted single-electron mesoscopic charge 
transport that was proposed in Ref. [9 . A number of 
the consequences of mechanically assisted "shuttle" trans- 
port have been investigated for normal, superconducting 
and magnetic NEM-SET devices (see, e.g., the reviews 
Ref. [TF and [TT). 

Other frequently studied nanoelectromechanical sys- 
tems are "one-dimensional" in the sense that they incor- 
porate a suspended nanowire as a (possibly) vibrating 
current carrying element. In this case there are many 
more relevant mechanical degrees of freedom than in the 
rigid "zero-dimensional" quantum- dot- type NEM devices. 
In the latter type of device it is sufficient to take the 
center-of-mass-motion of the movable element into ac- 
count, while in the former also fiexural and other types 
of wire vibrations, as w ell as int rinsic phonon modes, can 
couple to the electrons^^ The effect of the vari- 

ous vibration modes on electron transport can be investi- 
gated either through the driven resonator metho d, wh ere 
a down-mixing technique is used for detection^^^^, or 
through scanning tunneling spectroscopy methods, de- 
tecting phonon - assiste d channels for the inelastic tunnel- 
ing of electron^^^. 

A number of novel nanoelectromechanical phenomena 
become possible as a result of the elongation of the mov- 
able NEM-SET element. The present authors have, e.g., 
recently been involved in studies of the interplay in such 
devices between quantum coherence in both electronic 
and vibrational degrees of freedom. Other investigations 
have been related to non-equilibrium transport issues re- 
sulting from electromechanical instabilities that drive the 
system far away from equilibrium. This work will be re- 
viewed in the present paper, where the following topics 
will be discussed. 

1. Multi-mode shuttle structures. The nanoelectrome- 
chanics of nanowire-based ("one-dimensional") 
NEM-SET devices have to deal with the multimode 
fiexural vibrations of the suspended wire, and their 



effect on charge transport. The relevant theoreti- 
cal framework will be reviewed, and in particular 
we will show how a strong nonlinear coupling be- 
tween different modes leads to a "self organization" 
effect in multimode shuttle structures. 

2. Interplay between quantum coherent mechanical vi- 
brations and coherent electron transport. Quantum 
coherence is expected to play a significant role in 
the nanomechanics of sufficiently long suspended 
carbon nanowires, since the amplitude of their zero- 
point oscillations is comparatively large. The en- 
tanglement between quantum coherent electrons 
and quantum coherent mechanical vibrations, in- 
duced by an external magnetic field, will be shown 
to qualitatively modify the NEM-assisted electron 
transport in this case. This phenomenon suggests 
that measuring the magnetoresistance may be a 
way of testing quantum coherence in mechanical 
nanovibrations. 

3. Nanoelectromechanics of a superconducting weak 
nanowire-link. A current passing through a sus- 
pended wire of nanometer-sized cross-section in the 
presence of a transverse magnetic field, gives rise to 
a strong nanoelectromechanical coupling via both 
the Lorentz force and the electromotive force. We 
show that this coupling qualitatively changes the 
electrodynamics of a suspended nanowire serving 
as a superconducting weak link. The possibility to 
resonantly pump energy into nanowire vibrations 
by means of an AC Josephson current is explored. 

The paper is organized as follows. In Section [IT] we 
introduce the reader to the important concepts of this 
Review and in particular focus on the effects of charge 
concentration in nanoelectromechanical tunneling struc- 
tures. These ideas are expanded upon further in Sec- 
tion II A and II B| where the conditions for the single- and 
multimode shuttle instability respectively are derived. 
We conclude this part of the Review by demonstrating 
how self-selectivity between the vibrational modes of the 
considered suspended nanowire device can be achieved. 
Section III CI and IIIDI In Section IIIII we instead con- 
sider the effects of current variations coupled to an exter- 
nally applied magnetic field and analyze the subsequent 
electron-vibron interaction for the case of normal, Sec- 
tion |III A} and superconducting, Section [lIIB[ electrodes. 



II. ELECTROMECHANICAL COUPLING IN 
TUNNELING NANOSTRUCTURES WITH 
CHARGE CONCENTRATION 

Electron transport in tunneling nanostructures can be 
strongly affected by the accumulation of charge in small 
parts of the devices. The increase in electrostatic en- 
ergy corresponding to such charge accumulation can eas- 
ily exceed typical thermal energies and energies available 
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nanowire, the complete set of its fiexural vibrational 
modes should be taken into account. It will be shown 
later that, under certain conditions, this adds some new 
features to the ordinary picture of nanomechanical shut- 
tling in "point-like" movable islands. However, the con- 
ditions for the electromechanical instability to occur are 
not crucially affected by the multimode nature of fiexural 
vibrations, hence we will start the analysis of the shuttle 
instability considering only a single vibrational mode. 



A. Single mode instability 



FIG. 1: (Color online) Sketch of the system considered (a) 
and the equivalent circuit (b). Mechanical deviation of the 
wire away from the static configuration is described by the 
"displacement field" u{z, t). An STM tip is put over the point 
zo along the suspended nanowire axis and the electron tun- 
neling rate ri(u{zo)) — [Ri{u{zo))Ci{u{zo))]~^ between the 
STM tip and the nanowire depends on the spatial profile of 
the nanowire. The tunneling rate r2 = [R2C2]~'^ between the 
nanowire and the electrodes is, however, constant. Adapted 
with permission from 21 , L. M. Jonsson et ah, Nano Lett., 
5, 1165 (2005). © 2005, American Chemical Society 



from the voltage source, resulting in "Coulomb blockade" 
of electron tunneling^. The Coulomb forces associated 
with the uncompensated electric charge can furthermore, 
under the right conditions, be large enough to induce sig- 
nificant mechanical deformation of the movable parts of 
nanomechanical devices. The "electromechanical insta- 
bility" or "shuttle instability" of the NEM-SET device 
considered in Refs. [QilO'll is a remarkable example of 
this phenomenon. In this system electromechanical cou- 
pling is achieved due to the accumulation of charge on 
a movable metal island ("dot") and mechanical "shuttle" 
transport of single electrons can be achieved as a result 
of an electromechanical instability if a sufficiently large 
bias voltage is applied. 

In the following sections we will explore the possibility 
to utilize the fiexural vibrations of a suspended nanowire 
to transport electric charge. Such mechanically assisted 
transport can be attained if the electric charge injection 
is focused into the movable part of the device. This can 
be accomplished through the use of a scanning tunnel 
microscope (STM) which locally injects current into the 
vibrating suspended nanowire as shown in Fig. |l(a)[ 

An experimental realization of this device has been 
studied by LeRoy et al.^^ and details on the fabrication of 
the device can be found in Ref. [22 . In this work, the res- 
onator consists of a single-wall carbon nanotube (CNT) 
suspended over a trench such that a segment of the tube 
is free to move in response to external forces. Since the 
shuttle (in) stability analysis presented here does not cru- 
cially depend on the material of which the movable part 
is made, we will in the following consider a generic oscil- 
lating nanowire. 

In order to fully describe the mechanics of the 



The device shown in Fig. |l(a)] models the nanowire as a 
beam of length / with clamped ends. Let the undeformed 
wire extend along the z-axis, while its cross-section lies in 
the x?/-plane. If the wire is deformed from its static con- 
figuration by a bending force perpendicular to the z-axis, 
it will start to oscillate in the xz-plane. Such fiexural vi- 
brations of CNTs have been detected in devices similar 
to the one analyzed her d^^'^^'^^'^^ ' ^^ 1™. 

The deformation of the nanotube profile along the 
x-direction can be described by the displacement field 
ii(z,t), which, according to linear elasticity theor}^^^ 
obeys the following equation of motion. 
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Here, E is the Young's modulus of the nanowire, / is 
the cross-sectional area moment of inertia, p is the linear 
mass density, S is the cross-section area and / is the force 
per unit length acting on the nanowire. The boundary 
conditions for the doubly-clamped beam are ii(0,t) = 
u(l,t) = and du/dz{0,t) = du/dz{l,t) = 0. 

The displacement field and the force per unit length in 
Eq. ^ can be expanded in the complete set of normal 
modes {(Pj{z)} obtained by diagonalization of the her- 
mit ian operator d^/dz^ with the above boundary condi- 
tions. The result is 

^i^^t) = '^Xj{t)^j{^) (2a) 

j 

f{z,t) = ^fjiXi{t),X2it),...,t)^j{z), (2b) 

j 

where {xj{t)} and {/j(xi(t), X2(t), . . . , t)} are time- 
dependent amplitudes. Note that we write the ampli- 
tudes for the force as {/j(xi(t), X2(t), . . . , t)} instead of 
just {fj{t)} to stress that, in the system considered, 
the force at any time t depends also on the displace- 
ment of the wire, i.e. on the normal mode amplitudes 
{xj{t)}. Inserting the expressions ^ for u{z^t) and 
f{z^t) into the equation of motion ([T]) for the beam, one 
finds that the beam dynamics is equivalent to that of 
a set of coupled harmonic oscillators {xj} with frequen- 
cies {ujj}^ affected by forces {/j(xi, X2, . . . , t)}. In gen- 
eral, the unperturbed modal oscillation frequencies are 
given by coj = {cj /l'^){EI / pSY^'^ ^ where the coefficients 
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Cj = 22.4,61.7,120.9,199.9,298.6,... are obtained by 
solving the equation cos cosh ^cj = 1, see Ref. [25] 
for details. 

The normal mode eigenfunctions {^j{z)} are well ap- 
proximated by {sin(j7rz//)} and can be classified accord- 
ing to their symmetry properties; either "even" or "odd" 
under the spatial inversion operation z ^ —z with re- 
spect to the midpoint of the nanotube. The modes la- 
beled 1, 3, . . . , 2j + 1, . . . turn out to be even, while the 
modes 2, 4, . . . , 2j, . . . are odd. For the moment we as- 
sume that the STM tip is put over the midpoint of the 
nanowire {zq = 1/2). This assumption enables us to ne- 
glect the coupling between the STM and the odd modes. 
Furthermore, only the fundamental bending mode is 
taken into account, for which (pi{l/2) ~ 1 (if the normal 
modes {^j{z)} are normalized to 1). The displacement 
of the nanowire then reduces to the amplitude of the first 
mode (for the rest of this section indicated as x{t) instead 
of xi{t)). Finally, a DC bias voltage V is applied between 
the STM and the electrodes enabling electrons to tunnel 
from the STM tip to the nanowire and then from the 
nanowire to the electrodes. 

From the point of view of the charge transport the 
system can be modeled as two tunnel junctions in series. 
The junction between the STM tip and the nanowire, la- 
beled "1" in Fig. |l(b)[ is characterized by a capacitance 
Ci{x) and a resistance i^i(x), both of which depend on 
the defiection of the tube. This dependence is assumed to 
be of the form Ri{x) - g-^^+^^Z^ = i?oe"^/^, where A — 
the tunneling length — is a parameter used to character- 
ize the tunnel barrier. If the wire oscillates, its separation 
from the STM tip changes in time, hence so does the tun- 
neling probability. The interface between the nanowire 
and the metal electrodes is also assumed to be a tunnel 
junction (rather than an ohmic contact). It is labeled "2" 
in Fig. |l(b)| and is characterized by two constants, the 
capacitance C2 and the resistance R2. 

When an electron tunnels onto the nanowire, a certain 
time Tq is needed in order to redistribute the charge on it. 
This time can, however, always be assumed to be shorter 
than the other characteristic times in the system. These 



• Ti , the time needed for electrons to tunnel from the 
STM tip to the nanowire 

• T2, the time needed for electrons to tunnel from the 
nanowire to the electrodes 

• T, the period of nanowire oscillations. 

Assuming that Tq can be neglected allows us to treat the 
wire as a metal island with a well-defined excess charge 
q and an electrostatic potential When the wire is 
charged by electrons tunneling from the STM tip, it is 
affected by a capacitive force Fcap whose intensity de- 
pends on the voltage drop between the STM tip and the 
wire. 

The question is now under what circumstances a shut- 
tle instability can be expected to occur in the system, a 



problem that was initially addressed by Jonsson et al^. 
Here, we remind the reader that the shuttle instability 
requires that the static equilibrium state of the wire is 
unstable with respect to the formation of a new dynami- 
cal stationary state characterized by finite-amplitude os- 
cillations, so called limit cycle oscillations'^^ ^ around the 
static configuration. 

To answer the above posed question, let us consider 
the case where the nanowire is only slightly perturbed 
from its stationary configuration and let free to oscil- 
late. The only thing that can then happen, given that 
the bias voltage V is not too large, is that it performs 
some damped oscillations and returns to its initial con- 
figuration. In general a broad variety of microscopic 
mechanisms contribute to such damping of oscillations 
in nano-resonators. Some of them, e.g., thermoelastic 
damping^'^, are related to internal features of the res- 
onator while others, e.g., losses due to the clampin^^ 
and to air friction^^, depend on its interaction with the 
environment. For the moment we assume that in the sys- 
tem considered all dissipative effects can be taken into ac- 
count through a phenomenological viscous damping term 
—jx (the effects of other types of damping will be con- 
sidered in Section IIIDl). 



If the wire is instead perturbed further and allowed 
to perform larger mechanical vibrations the situation is 
somewhat altered. In this scenario the wire can be very 
close to the STM tip (-u < A), in which case we have that 
Ri <C i^2, i-e. Ti <C r2. Under these conditions tunneling 
from the STM tip is very likely and the electrostatic po- 
tential ^ of the wire is approximately the same as that 
of the STM tip. On the other hand, if the wire is in- 
stead far from the STM tip (li > A), the opposite is true 
and ^ approaches the potential of the electrodes. Now, if 
the wire oscillates with frequency cj, the parameters a;ri, 
Ljr2 have some finite values which implies that at any 
moment ^ has no time to adjust itself to the equilibrium 
value determined by the wire's position. As a result it 
also depends on the wire's position and velocity at earlier 
times. In other words, there is a correlation between the 
wire's velocity and the electrostatic force. This in turn 
implies that the force performs some work on the wire, 
(xFcap) 7^ 0, where (...) is the time average over one 
period. 

According to the ordinary shuttle theory, if the work 
done by the electrostatic force overcomes the work done 
by the dissipative forces, the electromechanical instability 
occurs. In order to check these qualitative considerations, 
the evolution of the system has been modeled in Ref. ^T] 
using the following equations of motion for the position 
of the wire's midpoint, x, and the excess charge, on it 



X -\- ^yx -\- ijx = 



Fo{C2V -qf _^ 

(1-f) 

• ^ fn f ^ g-(^) 



(3a) 



In Eq. ([3]), Cy.{x) = Ci{x) + C2 is the total capacitance. 
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FIG. 2: Threshold dissipation as a function of the junction 
resistances {Cf. Fig. |l(b)] ). The largest values of 1/Qthr are 
obtained when the two junction resistances are of the same 
order of magnitude. Reprinted with permission from [21 , 
L. M. Jonsson et al, Nano Lett, 5, 1165 (2005). © 2005, 
American Chemical Society. 



C-{x) = Ci{x) -C2,G^ = l/Ri{x) + I/R2 is the total 
conductance and G-{x) = 1/Ri{x) — l/i^2- Due to the 
short distance between STM tip and nanowire, surface 
forces might also be significant. However, their relevance 
compared to the elasto-mechanical forces can always be 
reduced by a suitable choice of length and thickness of the 
wir^. The excess charge on the tube, is treated as a 
continuous variable, an approximation which is justified 
if the tunneling times ri, r2 are much smaller than the 
period of oscillation and the temperature is not too low. 

In order to perform a stability analysis of ([s]), one can 
linearize the two equations around the stationary solu- 
tion (xo, ^o)- From this analysis it is found that if the 
dissipation coefficient exceeds a certain threshold value, 
7 > 7thr for a fixed bias voltage V, the stationary so- 
lution is a stable fixed point. Therefore, any trajectory 
that starts close to this solution will ultimately fall into 
this point. However, if instead the dissipation coefficient 
is smaller than the threshold value, 7 < 7thr, any small 
initial deviation from (xq, ^0) grows in time and the sys- 
tem manifests a "shuttle-like" instability. The conditions 
for this instability can be reformulated as 1/Q < 1/Qthr, 
where Q is the "quality factor" of the oscillator. This is an 
experimentally accessible parameter that expresses the 
robustness of the oscillations against all possible sources 
of dissipation, hence it can be measured independently 
on the specific model used to describe dissipation in the 
system. For the case considered here, the Q-factor can 
be expressed as Q = I/7. Values of Q between 10 and 
10^ have been reported for carbon nanotube-based nano- 
resonators^^ and of the order of 10^ for SiC nanowire- 
based oscillators^^. 

In Fig. [2] we show the threshold dissipation coefficient 
7thr as a function of the parameters Rq, R2 (for a fixed 
bias voltage V). From the shape of 7thr we can extract 
several bits of information about the physical conditions 
under which the shuttle instability occurs. First of all, a 



negative or zero value for 7thr means that no instability 
can be established. As explained above, some retarda- 
tion effect (present for any finite value of uoti and UOT2) is 
necessary in order for the instability to occur, as the net 
work of the electrostatic force over one period of osciha- 
tion would otherwise be zero. The smaller the size of the 
retardation the more difficult it is to make the system 
unstable. This is also confirmed in Fig. |2j where 7thr ap- 
proaches zero for small values of the resistances Rq and 
R2 (implying uoti^ UOT2 0). On the other hand, it is 
clear that if the tunneling of electrons becomes rare due 
to large junction resistances, the nanowire is too weakly 
affected by the electrostatic force to make the instability 
occur. This observation is consistent with the behavior 
of 7thr in Fig. [2j where 7thr ^ in the limit of large 
resistances. Thus, we learn from studying Fig. [2] that 
the optimal condition to achieve the shuttle instability is 
to have the junction resistances very close to each other, 
Rq ^ R2' This implies uoti ~ UOT2 ~ 1, i.e. quite substan- 
tial retardation effect. Having Q-factors between 10 and 
1000 and a bias voltage of the order of 1 V, it seems pos- 
sible to switch on the instability for a large set of parame- 
ters Rq^ R2. In the experiment performed by LeRoy et al. 
the junction resistances are quite different {Rq ^ i?2^ 
which makes the electromechanical instability unlikely. 
For very asymmetric junctions, the electrostatic poten- 
tial of the wire at any time is mainly determined by the 
strength of the coupling to the lowest resistance. 

The physical condition for the occurrence of the shut- 
tle instability can also be expressed in terms of the bias 
voltage, instead of the dissipation coefficient, 7. Natu- 
rally this is an experimentally more meaningful quantity, 
as the amount of dissipation in the system is difficult 
to evaluate through only a single parameter, while the 
bias voltage can be externally controlled. The evolution 
of the amplitude of oscillation as a function of the ap- 
plied bias voltage (for a fixed 7, i.e. fixed quality factor 
Q) is showed in Fig. 3(a)[ 
ence of a threshold vo' 



From this image, the pres- 
tage above which the instability is 
established is clearly visible. The corresponding I — V 
and dl/dV — V characteristics have also been calculated 
according to the numerical solution of the equations of 
motion and are shown in Fig. |3(b)[ where the insets refer 
to the different ways in which the transition to the new 
equilibrium state may occur. This transition depends on 
the bias voltage and the parameters of the system, and 
can be classified as either smooth ("soft" excitation) or 
step-like ("hard" excitation)^^. These plots provide indi- 
cations of what is expected to be found in an experimen- 
tal investigation of the shuttle instability in suspended 
nano-resonators . 



B. Multimode shuttling of single electrons 

We now generalize the analysis of the system to include 
the full multimode description of its mechanical degrees 
of freedom. The formal similarity between the standard 
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FIG. 3: (a) Amplitude of nanowire oscillations as a function 
of the bias voltage. The solid line shows a "soft" develop- 
ment of the instability, while the dotted line corresponds to 
a "hard" transition, (b) I — V curves for the same "soft" 
(left) and "hard" (right) transitions. Insets show correspond- 
ing dl/dV — V curves. Reprinted with permission from [21 , 
L. M. Jonsson et al, Nano Lett, 5, 1165 (2005). © 2005, 
American Chemical Society. 



gies of the system (i.e. eV and kBT). There is a "contin- 
uum" regime, defined by eV, ksT ^ A, where quantiza- 
tion is not important, and a "discrete" regime, defined by 
eV, ksT <C A, where the electronic states on the tube are 
so far apart in energy that only one of them falls within 
the window set by the bias voltage. Under the latter 
conditions, only one state is involved in the electronic 
charge transport. As compared to the previous section, 
we now describe the tunneling junctions through the tun- 
neling rates Ti and where the subscripts "1" and "2" 
refer respectively to the junction between the STM tip 
and the wire and to the junction between the wire and 
the electrodes. The right tunneling rate, is assumed 
to be constant, while the left rate, Ti, depends on the 
tube displacement. What differs the "continuum" and 
"single-level" regimes is that the tunneling rate is voltage- 
dependent for the former, whereas this is not the case for 
the latter. Furthermore, we assume that the tempera- 
ture is low enough that "backward" tunneling processes 
(from the right to the left junction) are not allowed and 
charging of the tube with more than one electron is pre- 
vented by Coulomb blockade {ksT < e^/2C, where C is 
the capacitance of the nanowire). 

In order to study the motion of the system it is conve- 
nient to give a quantitative definition of the electrome- 
chanical coupling of the different mechanical modes. Let 
us consider the wire initially at rest in its stationary con- 
figuration, with no extra charge on it and let L be the 
distance from the STM tip. As soon as the bias voltage 
is switched on and an electron tunnels onto the wire the 
electrostatic field produced by the STM exerts an attrac- 
tive force on it. The nanowire moves from its stationary 
configuration with the corresponding displacement, Suq, 
being given by the equation of motion derived from linear 
elasticity theory, Eq. ([T]), 

^^o = ^E^- (4) 



shuttle model and the approach used here allows us to in- 
clude also Coulomb blockade phenomena, electronic level 
quantization and quantum description of the dynamics. 

In order to make the system sketched in Fig. |l(a)| 
equivalent to a single electron transistor device we will 
now consider a different bias voltage for the two elec- 
trodes. The left electrode is assumed to be at potential 
— y/2, while the right is capacitively connected such that 
it forms a gate electrode with potential Vg. No tunneling 
is allowed between the nanowire and the right electrode 
as in a standard SET device. The analysis of this model 
was first presented by Jonsson et alW^. 

Considering that the nanowire has a finite length, the 
set of single-particle states available for electrons tunnel- 
ing on it is quantized. We assume that these states are 
equally spaced in energy by the amount A. Under these 
conditions, two regimes can be distinguished according to 
the size of A compared to the other characteristic ener- 



Here, the sum is extended to all modes whose shape 
profile is antisymmetric with respect to the wire mid- 
point (j = 1, 3, 5, . . .). For larger absolute values of 6uo 
and smaller separations between the nanowire and the 
STM tip the probability of tunneling increases apprecia- 
bly. The characteristic length that provides a natural 
reference for all the other lengths in the system is the 
tunneling length A. Therefore, the ratio between the dis- 
placement due to the electrostatic force and the tunneling 
length Suq/X = e, provides an estimate of how much the 
motion of the wire is affected by the tunneling of charge. 
The ratio between each term in the sum in Q and A can 
thus be interpreted as the electromechanical coupling for 
each mode. The electromechanical coupling strength s 
is the same parameter that defines the ratio between the 
electrostatic force and the elastic force in the theory of 
the ordinary shuttl^. 

Since the shape profiles of the normal modes {^j{z)} 
are oscillating functions of the coordinate z whose wave- 
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length decreases with increasing mode index j, the cou- 
phng of the higher modes to the STM tip is weaker 
than to the lower modes. Therefore, a valid approxi- 
mation is to take into account only a limited number 
K of flexural modes. In the quantum description de- 
veloped in Ref. [32 , the electromechanical coupling for 
the j— th mode is defined as eljS /hwjXj = dj/Xj^ where 
Ij = ^/h/m^Jj and Xj are respectively the zero-point am- 
plitude and the effective tunneling length for the j— th 
mode. The electromechanical coupling can be defined 
and evaluated also for other types of mechanical modes, 
as shown in several works, Refs. |15|18|19|23j . 

The dynamics of this system can be described through 
a generalized master equation for the two reduced den- 
sity matrices that describe the neutral, po, state of the 
nanowire and the state with one extra electron, pi (higher 
charge states are not allowed because of the Coulomb 
blockade), 

^ = -i[Hosc,Po] -^^Rpi - ^{Tl{^,V)}pi 

~ 2 Z^(^[^^' ^o}] + [^^' [^^' ^o]]) '^^^) 

n 

^ = -i[Hosc{^ - d),/9l] + ^/^L{±,V)po^/TL{±,V) 

-r^Pl - ^'^{i[x„,{nn,Pl}] + [Xn,[Xn,Pl]]) ■ (5b) 

n 

In Eq. ([5|, Hose is the quantum mechanical Hamiltonian 
of the vibrational modes of the wire that, as in the clas- 
sical treatment, form a set of non- interacting harmonic 
oscillators, 

n ^ ^ 

where Pn = ^n/^i and {Tr^} are the momenta conju- 
gated to the mode amplitudes {f n}- The charging effects 
and the electromechanical coupling between the nanowire 
and the STM tip enter into the description of the system 
through the term 

^tube = ^ E{Vg)clck + -^n^ - hUifl ^ XjEj , (7) 
k j 

with E(Vg) being the single electron energies, [ck] 
are creation [annihilation] operators of single electron 
states labeled by the momentum U is the electro- 
static energy of the nanowire and Sj is the electromechan- 
ical coupling for the j— th mode. Dissipative effects are 
here assumed to be relatively weak and are characterized 
by a frequency-independent dissipation rate coefficient 7 
("Ohmic damping"). 

From Eq. ([5|, average displacements, momenta and 
probabilities can be found from the reduced density ma- 
trices and their respective equation of motion can be de- 
rived. These equations can be linearized around the sta- 
tionary solution and the behavior of the solution can be 
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FIG. 4: (Color online) Threshold dissipation for the first 
(7thr,i) and third (7thr,3) transverse mode as a function of Fq. 
Going from top to bottom the electromechanical coupling is 
increased from the weak to the strong regime, di/X < 1 (a), 
di/X ^ 1 (b) and di/X ^ 1 (c). When the electromechani- 
cal coupling is large the interaction between different modes 
is no longer negligible, which introduces a qualitatively new 
feature compared to the weak electromechanical regime. As 
can be seen, some values of Fo allow for mode n to be ex- 
cited and mode m to be suppressed even though n > m, i.e. 
7thr,n > 7thr,m- Reprinted with permission from [32], L. M. 
Jonsson et al, New J. Phys, 9, 90 (2007). © (2007), Deutsche 
Physikalische Gesellschaft. 
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investigated in order to check the onset of the instabihty. 
In the weak electromechanical coupling regime this anal- 
ysis can be carried out analytically and an expression for 
the threshold dissipation can be found. The results ob- 
tained indicate that the different modes can be treated 
independently from each other, as is shown in Fig. |4(a)| 
The strong electromechanical coupling regime must how- 
ever by analyzed numerically. The consequences of the 
transition from the weak to the strong coupling are shown 
in Fig. |4(b)| and |4(c)] where the threshold dissipations for 
mode 1 and 3 are plotted as functions of the coefficient 
of the left tunneling rate Fq. These curves demonstrate 
the possibility, in the strong coupling regime, to excite 
a certain mode m without also making all the modes 
n < m unstable, something that is unavoidable in the 
weak regime. The parameter Fq is experimentally con- 
trollable as it is determined by the distance between the 
STM tip and the nanowire in its stationary configura- 
tion. This distance can be measured and controlled with 
sufficient accuracy to make realistic devices. 



A natural issue that needs to be clarified is what hap- 
pens to the system once the instability develops. From 
the theory of the ordinary shuttle it is known that the 
amplitude of the oscillations should increase and then 
saturate to a value determined by the parameters of the 
device, the system reaches a limit cycle. In order to check 
the existence of the limit cycle for the system consid- 
ered here, linearization of the equations of motion is not 
enough and the full nonlinearity due to the exponential 
form of the left tunneling rate must be taken into ac- 
count. This issue has been investigated in Ref. [32] for 
the case of weak coupling and a single unstable mode. 
The existence of the limit cycle can be proved using the 
Wigner representation of the density operator, following 
the approach developed in Ref. [33]. The advantage of 
this approach is that it makes the crossover between the 
tunneling and shuttling regimes more transparent, as dis- 
cussed in Ref. [3T. From this analysis it is found that 
until the electromechanical instability develops, the sys- 
tem behaves basically like a series of two tunneling junc- 
tions, hence the only expected contribution to the current 
comes from tunneling. However, once the instability sets 
in and the system reaches its limit cycle (characterized 
by steady amplitude oscillations) the current is drasti- 
cally modified. In Fig. [5] the current is plotted as a func- 
tion of the electrostatic energy eV for the case of a single 
mode and "hard" instability transition . As can be seen 
in this image the characteristics of the "hard" transition 
are clearly visible in the sudden increase in the ampli- 
tude (and hence current) when the bias voltage exceeds 
the threshold value. Another distinctive feature of this 
process (that is absent for the case of "soft" transition) 
is hysteresis. Sweeping the voltage down from a value 
above the threshold, the current does not correspond to 
the values found before the onset of the instability. 



0.14 



0.1 



0.06 
70 



- - - Decreasing voltage 
^Increasing voltage 


* 

* 

* 

i ■ 

1 









75 



80 

eVifiG)) 



85 



90 



FIG. 5: (Color online) / — ^/-characteristics for the device of 
Fig. 1(a) in the weak coupling regime. For small voltages, the 
current is constant and depends only on the tunneling through 
the double static junction / = eVoV n/Vt. Above the thresh- 
old voltage the current also depends on the vibration ampli- 
tude. Here the electromechanical instability occurs through 
the "hard" transition, characterized by the displayed hystere- 
sis in the / — V-curve. Reprinted with permission from [32], 
L. M. Jonsson et al, New J. Phys, 9, 90 (2007). © (2007), 
Deutsche Physikalische Gesellschaft. 



Mult ist ability and self-organization of 
multimode shuttle vibrations 



The consequences of the electromechanical instability 
on the whole set of flexural modes (not only on the low- 
est frequency mode) has been investigated by Jonsson et 
al.^^ ^ following a classical approach. Starting from linear 
elasticity theory, the equations of motion for the ampli- 
tudes of the transverse normal modes were found to be 
(the STM tip is still supposed to be at //2 so only even 
modes with respect to the midpoint of the tube are in- 
cluded), 



■i^lxn = q£/m. 



(8) 



Here, the force term that acts on the wire depends on the 
average excess charge, on it and the effective electric 
field E produced by the applied bias voltage between the 
STM tip and the electrodes. The dissipation is modeled 
by a linear viscous term —jXn and, since the damping 
coefficient 7 is assumed to be constant, every mode is 
damped in the same way. 

An important parameter for the system under consid- 
eration is the ratio between the typical mechanical vibra- 
tion frequency u and the characteristic rate of electron 
tunneling F. In the following we assume that tunneling 
events occur frequently compared to the period of oscil- 
lations, i.e. a; <C F. Even though this is far from being 
the optimal condition for the instability to occur, as dis- 
cussed in Section [ITAI it allows us to write down a simple 
kinetic equation for the time evolution of the probability 
of having one extra charge on the wire (see below). This 
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choice will result in more experimentally demanding con- 
ditions to achieve the instability (i.e. higher Q-factors) 
with respect to those previously discussed in Section 
A more general approach will be followed in Section 

As was done in Ref. [32 , we here assume that the tem- 
perature, the bias voltage and electrostatic charging en- 
ergy of the nanowire are such that backward tunneling 
processes are forbidden and no more than one extra elec- 
tron is allowed on the nanowire at any given time. Let 
us indicate with pi{t) and po{t) respectively the proba- 
bility to find one and zero electron on the nanowire at 
time t. Then, the average excess charge at time t can 
be expressed as q{t) = epi{t). The variation of pi{t) in 
time is given by the difference between the probability 
of tunneling from the STM tip to the nanowire when the 
latter is neutral and the probability of tunneling from the 
nanowire to the electrodes when there is already one ex- 
tra electron on it. We also have that the probabilities for 
the allowed charge states must sum up to 1 (pi +po = 1, 
higher charge states are prevented by the Coulomb block- 
ade). Exploiting this fact, we can write down the rate 
equation for p{t) = pi{t) in the form. 



p{t) = -r{uo)p{t)^ri{uo), 



(9) 



where T{uo) = Ti{uo) + r2 with uq = Hn^n^nl^o) ~ 
S^x^, while Ti{x) = Fq exp(— x/A) is the rate of electron 
tunneling across the STM-nanowire junction, so that the 
typical rate of electron tunneling is F = Fq + The 
solution of Eq. ^ can be written as a series expansion 
in the parameter coi/T which is assumed small. 

Equations ^ and ^ describe the coupled nanoelec- 
tromechanical dynamics of the nanowire-based NEM- 
SET device. Since strong bias voltages are not allowed 
if the excess charge on the wire must be limited to one 
electron, it is reasonable to consider the limit of weak 
electromechanical coupling. In this regime the displace- 
ment of the nanowire due to a single excess charge is small 
on the scale of the tunneling length. Then, as shown in 
Ref. [32 , the onset of the shuttle instability occurs inde- 
pendently for each vibrational mode and if c<; <C F the 
instability is "soft' ^^. 

In order to analyze the evolution of the system after 
the instability occurs, it is not sufficient to linearize the 
equations of motion ([8| and ([9| around the stationary 
configuration as the linear force term in the right hand 
side of these equations simply tells if the instability can 
develop or not. Therefore, it is necessary to take into 
account at least some nonlinear terms. Since within the 
regime of weak electromechanical coupling the amplitude 
at which the nanowire oscillates is expected to be small 
compared to the tunneling length, it is possible to expand 
the tunneling rate Fi(iio) for small uq/X and keep con- 
tribution up to the third order. Another approximation 
that is physically motivated in the weak coupling regime 
is the use of the Ansatz Xn{t) = XAn{t) sm{ujrit + Xn{t))^ 
where the amplitudes An{t) and phases Xn{t) are slowly 
varying functions on the time scale defined by the pe- 
riod of oscillations. This difference in characteristic time 




FIG. 6: (Color online) Stationary points when two modes 
(n and m) are unstable. Two attractors, indicated by (o), 
corresponding to a finite amplitude of one mode while the 
other mode is suppressed are shown. The stationary point 
marked with (x) is a repellor and the point indicated by (*) 
is a saddle point. The thick lines are separatrices that trajec- 
tories cannot cross. The separatrix An = Am ensures that if 
An{0) > Am(0), this inequality hold for all times t. Reprinted 
with permission from ^5]? L. M. Jonsson et ah, New J. Phys., 
9, 90 (2007). © (2007), Deutsche Physikalische Gesellschaft. 



scales allows us to replace the Ansatz into Eq. ^ and 
take the time- average over the period of the first, i.e. 
longest, mode. In this time interval the amplitudes and 
phases for the different modes can be considered constant 
and be replaced by their average values. Furthermore, 
the time variation of the phases is assumed to be negligi- 
ble with respect to the eigenfrequencies of the nanowire, 
Xn cjn, which means that only the amplitudes An are 
taken into account in the description of the system. The 
resulting equation of motion for the averaged amplitudes 
are. 



An — OLnAn I ^ ^ ^ A^ 

ra^n 



(10) 



where 5n and an are defined from the following combi- 
nations of parameters 



128Ar 



l + O 



r2 



(lib) 



To first order in cj^/F the expressions for 5n and do 
not depend on the mode index n, as the product uo'^dn 
is independent on n, hence one can replace 5n ^ 5 and 



The behavior of the solution of Eq. ( 10 ) can be explic 



itly visualized for the case of two modes n, m (the gener- 
alization to more modes is straightforward). The corre- 
sponding stationary points of the two nonlinear coupled 
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equations can be found analytically and their stability 
can be determined through the evaluation of the Jaco- 
bian matri:jP^. As the dynamical behavior of the system 
depends on the sign of the parameter 5 we identify, from 
Eq. (p^Ol), that if J < the only stationary point is the 



origin, corresponding to the absence of any oscillation. 
For this set of parameters the nanowire is at rest in some 
static configuration determined only by the constant tun- 
neling rates Fq and F2. On the other hand, if 5 > 
the origin becomes unstable and three more stationary 
points appear: a saddle point at (^/3,^/3) and two sta- 
ble points at (aA, 0) and (0, vA). These two new stable 
points represent oscillating states with finite amplitude 
and frequency cOn or cOrn respectively. Which of these 
two will be reached by the system depends on the initial 
conditions as shown in Fig. [6j 

The conditions 5 > defines the onset of the shuttle 
instability. From the solution to Eq. ([lO| with (5 = it 
is possible to find the expression for the threshold elec- 
tric field above which the instability starts to develop, 
£c = 4TX^m/e. Further analysis of Eq. (10) indicates 



that once the instability for a certain number of vibra- 
tional modes is established, the system evolves in such 
a way that only one of the unstable modes reaches the 
new stationary state, characterized by steady amplitude 
oscillations, i.e. the limit cycle. It is also found that the 
selection of the surviving mode is determined by the ini- 
tial conditions as the mode which initially has the largest 
displacement from the origin (that is from the static equi- 
librium state) maintains its separation from the other 
modes and evolves into the limit cycle. This can be un- 
derstood as trajectories in the amplitude space cannot 
cross, a result that can be analytically generalized to an 
arbitrary number of modes by studying the asymptotic 
behavior of the solutions of Eq. (p!Q|^^. 

It is worth to remark that the symmetry between the 
modes that characterizes Fig. [6] is actually broken if 
higher order (in oj^/F and com/^) corrections to S be- 
come relevant or if the dissipation affects each mode in a 
different way. Then each mode will have its own Sn and, 
in general, the selection of the surviving mode may not 
only depend on the size of the initial displacements from 
the static equilibrium point. 

In order to check the results obtained by the analysis 
of Eq. (10), one can also compare them to the numer- 
ical solution of the "full" equations of motion, Eq. ([8|. 
The comparison for a given choice of initial conditions 
is shown in Fig. |7(a)[ Since the frequencies of the dif- 
ferent modes are not commensurable, i.e. they are not 
integer multiples of the fundamental mode, the initial 
motion of the nanowire is not characterized by a sharply 
defined periodicity as can be seen in the lower left panel 
of Fig. |7(a)[ However, when the system reaches the final 
stationary state and only the mode n = 5 (which initially 
had the largest deviation away from static equilibrium) 
has a non-zero amplitude, its oscillations are clearly peri- 
odic with frequency co^. Similar figures with modes n = 1 
or 3 having finite amplitudes can be obtained by changing 
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FIG. 7: (Color online) Numerical solution of (|8| and ([9| for 

the nanotube vibration amplitude as a function of time when 
three modes (n = 1,3,5) are unstable (S > 0). (a) Weak 
electromechanical regime. Comparison with the approximate 
result, Eq. (10), is shown as dashed curves. The lower left 



panel shows the quasi-periodic oscillation of the nanotube cen- 
ter position just after the onset of the instability, while the 
lower right panel shows the regular vibrations that appear af- 
ter all but the n=5 mode amplitudes have been suppressed 
(see text), (b) Strong coupling regime. The large amplitudes 
make an approximate analysis based on Eq. (10) invalid, but 
the phenomenon of selective excitation persists. Reprinted 
with permission from [35 , L. M. Jonsson et al, New J. Phys., 
9, 90 (2007). © (2007), Deutsche Physikalische Gesellschaft. 



the initial conditions. 

Note that the numerical analysis of Eq. ([s]) is not nec- 
essarily limited by the requirement that the oscillation 
amplitude be kept small with respect to the tunneling 
length. Fig. |7(b)[ for instance, suggests that the se- 
lective evolution promoted by the simultaneous instabil- 
ity of many vibrational modes can characterize also the 
"large amplitudes" regime. 
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D. Geometrical scanning of the flexural modes by 
STM tip displacement 

The theoretical study described in the previous section 
suggests that the development of the shuttle instability 
in an extended object, such as a suspended nanowire, 
induces a selective evolution of its many mechanical de- 
grees of freedom. It should be stressed however, that the 
results presented in Ref. [35 provide more a theoretical 
demonstration of this selectivity rather than an experi- 
mental procedure to control it. The prime reason for this 
is that the analysis involves the separate choice of initial 
conditions for each vibrational mode, an operation that 
cannot be performed experimentally. However, the de- 
vice considered can still be used to accurately detect the 
instability of the survivor mode, as the shuttling current 
is proportional to the mechanical frequency once the sys- 
tem reaches the stable limit cycle^. Thus, by measuring 
the current, we are able to tell which vibrational mode 
has become unstable. 

Besides this, the question of finding a practical proce- 
dure to select single vibrational modes through the shut- 
tle instability naturally arises. A possible way to do this 
was suggested in Ref. [36 , by generalizing some of the 
features of the system considered in Refs. j 21|32|35j . In 
this work the tip of the STM is specified in a generic po- 
sition zq along the nanowire axis and is no longer fixed 
above the wire's midpoint. Under these circumstances 
the coupling e is affected by the shape profiles of the 
eigenmodes at point zq (see Eq. Q). 

Ref. [36] also considers the effects of dissipation in the 
studied device to greater detail than what has been pre- 
sented in the previous sections. Damping of oscillations 
in solids can be caused by many different microscopic 
mechanisms and in general it is impossible to formulate a 
theory that is able to describe all of them from first prin- 
ciples with arbitrary accuracy. In most cases it is even 
difficult to accurately determine which is the dominant 
mechanism for dissipation. For example, certain dissi- 
pative processes that can be neglected in bulk materials 
can be enhanced in nano-devices due to the increased 
surface-to-volume raticP^. To account for these consid- 
erations Ref. 2^ analyzes not only the usual "viscous" 
damping term —^ii but also considers the Zener theory 
of dissipation for the standard linear solicP^. 

The Zener model provides a simple approach that goes 
a step beyond the ordinary elasticity theory. Accord- 
ing to Hooke's law, the stress and strain fields, i.e. the 
macroscopic variables that characterize the mechanical 
state of a solid, are connected by a simple proportional- 
ity relation, such that any change in the former reflects 
itself instantaneously as a change in the latter and vice 
versa. In the Zener theory, it is instead assumed that 
these quantities need a finite relaxation time, and 

respectively, to reach their equilibrium value. This 
is a consequence of the internal processes that dissipate 
energy during any modification of the mechanical state 
of the solid. For the case of periodic variations of the 



stress and strain with frequency cj, the following nonlin- 
ear frequency-dependent expression for the Q-factor can 
be derived from the Zener theory. 



1 1 + (c^f)2 



(12) 



Here, f = s/r^ and A = Eu — Er is the difference be- 
tween the "unrelaxed" and "relaxed" values of the Young's 
modulus^^, in effect a measure of the degree of non- 
elasticity of the body. 

In Ref. [36j a cla ssical description of the system de- 
picted in Fig. |l(a)| was formulated where the tunneling 
process that charges the nanowire and makes it sensi- 
tive to the electrostatic force is a stochastic process. As 
such, the normal mode amplitudes {xn} and conjugated 
momenta {tt^} are also stochastic variables and we de- 
fine probability densities for them when the nanowire 
is charged tti, X2, 7r2, . . . , t) and when it is neutral 

Po(^i, TTi, X2, 7r2, . . . , t). The evolution of these objects is 
determined by a generalized Boltzmann equation where 
the "collisional integral" is replaced by tunneling terms 
following the approach used in Ref. ^39j , 



dPi 



dPp 
dt 
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(13) 



(14) 



The behavior of the system is described by the dynamical 
variables obtained by averaging the mode amplitudes and 
conjugated momenta over Pq and Pi: (...)+ = (.. .)o + 

The possibility of having a shuttle instability in the 
weak electromechanical coupling regime is investigated 
by linearizing the equations of motion for the averaged 
dynamical variables around the static equilibrium po- 
sition where the fundamental solutions have the form 
~ e^^^^^^ . The condition for instability of the 
static configuration of the j-th mode can then be ex- 
pressed as, ^e[(^j] > 0. For a given quality factor Qj 
(that is, for a given amount of dissipation), the expo- 
nent 5j for each mode can be plotted as a function of 
the two experimentally accessible parameters: the tun- 
neling rate Fq between the STM tip and the nanowire in 
the static configuration and the position of the STM tip 
along the wire axis. The result of this analysis is plotted 
for t he case of viscous dissipation {Qj ex uoj) 
8(b) I for the case of internal dissipation. The 
two plots. Fig. |8(a)] and |8(b)[ do not show any qualita- 
tive difference. In both cases one can distinguish sets of 
parameters (Fq, z) for which only one mode is unstable 
and sets for which two or more modes are unstable. This 
fact suggests a general way to express the condition for 



in Fig. 8(a) 
and in Fig. 
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FIG. 8: (Color online) (a) Regions of instability for dissipa- 
tion modeled by a viscous term —^ii. (b) Regions of insta- 
bility for internal dissipation given by the Zener model. The 
solid line and the dots that define the threshold dissipation 
curves are obtained respectively from the analytic and the nu- 
merical solution of the linearized equations of motion. In both 
plots the filled areas define the values of z and To for which 
only a single mode is unstable. Reprinted with permission 
from |36|, F. Santandrea 



shuttle instability for other types of mechanical degrees 
of freedom and dissipative processes, ^e[(5j] > 0, where 
the real part of 5j is given by, 

_ ^3 , ToF^ eS ^]{z,) 
' 2Qj 2F, mAF,c.|+Fr 

III. MAGNETIC FIELD INDUCED NEM 
COUPLING 

In the previous sections we have shown that NEM- 
SET shuttle structures can be used to significantly 
alter the electronic and mechanical characteristics of 



the discussed system. In what follows, we will consider 
further the case of suspended nanowire structures and 
show that these can also be made to mechanically 
oscillate through the introduction of external magnetic 
fields. Due to the effective induced electron-vibron 
coupling of these devices novel physical phenomena are 
predicted as presented below. 

Although suspended nanowires can be viewed as a par- 
ticular realization of a NEM-SET shuttle structure, spe- 
cific nanoelectromechanical operation is expected as a re- 
sult of strong elongation of the movable part of these de- 
vices. This comes about as the very large aspect ratios 
in these systems affect both the mechanics of the flexu- 
ral vibrations and the electrodynamics of the electronic 
current flow. Elongation of a suspended wire, or alter- 
natively diminishing its cross-section, will make the wire 
more flexible and hence more sensitive to external me- 
chanical perturbations. This will in turn make quantum 
effects more pronounced due to the larger amplitude of 
the wire's zero point flexural vibrations. Utilizing this, 
new types of NEM coupling can be achieved through the 
induced electron-vibron interaction caused by, e.g., ex- 
ternally applied magnetic fields. This coupling specifies 
the Lorentz force acting on the wire, and for case of me- 
chanically oscillating wires, also the electromotive force 
on the electrons which counteracts the motion causing 
the vibrations. Since the flexibility of the wire increases 
with its aspect ratio, and the Lorenz force is directly pro- 
portional to the length of the wire, significant inductive 
NEM coupling can be achieved as a result of strong elec- 
trical current concentrations in wires of nanosized cross- 
section. 

In comparison to the nanoelectromechanics considered 
in the previous parts of this paper — where the tuning of 
the NEM performance was achieved through the coupling 
of an external electric field to the local charge concentra- 
tion — we here show that similar results can be achieved 
if external magnetic fields coupled to the electronic cur- 
rent are instead used. The effects of this type of NEM 
coupling on a suspended nanowire of length L and cross- 
section 5, carrying a current / in a transverse magnetic 
field H is described through the induced deflection of the 
wire, Xd{H). From standard elasticity theory one finds 
that Xd{H) oc Hj/Sa{L) where a{L) = E/L^ with j the 
current density and E the Young's modulus of the wire. 
For suspended nanowires of length L ~ 1 /im we estimate 
that such deflections can be as large as Xd{H) ~ 0.1- 
1 nm, which can crucially affect the electronic tunneling 
through such mesoscopic NEM structures. Quantum co- 
herence in the flexural vibrations of the wire may, as such, 
affect the electronic transport in non-trivial ways if the 
area Lxq {xq is the zero point amplitude of oscillation of 
the wire) available for penetration of the external mag- 
netic field is comparable to the flux quantum. For the re- 
alistic case of a vibrating carbon nanotube, such quantum 
coherent nanoelectromechanics can be achieved for mag- 
netic fields of the order of a few tens of tesla, correspond- 
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FIG. 9: Schematic diagram of the system considered in Sec- 
tion [mA] A transverse magnetic field, H, is appHed to a 
suspended 1-dimensional model nanowire of length L. If the 
wire is biased by a voltage V, it carries a current and the wire 
oscillates in response to the induced Lorentz force. Quantum 
fluctuations in the wire's bending modes make the electrons 
propagate along an effectively 2-dimensional wire. The mag- 
netic field-induced reduction of electron propagation in the 
elastic channel can be interpreted as an effect of destructive 
Aharonov-Bohm-type quantum interference between different 
paths of the tunneling electrons. Together with a blockade 
of some inelastic channels due to Pauli-principle restrictions 
(see text and the caption to Fig. [lO| this leads to a finite 
magnetoresistance of the wire. Amplitude shown is greatly 
exaggerated. Reprinted with permission from [40], G. Sonne. 



ing to present state of the art experimental achievements. 
In the following sections two examples of the above dis- 
cussed inductive nanoelectromechanics will be presented. 
The new physics of entanglement of resonantly tunnel- 
ing electrons with nanowire quantum vibrations will be 
presented in the Section [Tll A[ whereas the possibility to 
pump nanovibrations by a supercurrent flowing through 



a suspended nanowire is discussed in Section IIIB 



A. 



Quantum mechanically induced electronic 
Aharonov-Bohm interference 



Quantum fluctuation of the nanowire bending modes 
is a feature which must be taken into account while con- 
sidering the effects of quantum coherence in the vibration 
dynamics of a suspended nanowire based tunneling device 
(see Fig. [9|. This implies that geometrical constraint for 
the electrons, set by wire geometry, will no longer localize 
them to the 1-dimensional conducting wire, but will be 
also imply certain delocalization in the direction trans- 
verse to the wire axis. Such change of the dimension- 
ality opens up the possibility for quantum interference 
effects in the electronic tunneling over the leads if an ex- 
ternal magnetic field, is applied. As a result, a finite 
magnetoconductance of the ID vibrating wire occurs as 
a manifestation of nanomechanical quantum coherence 
eflFects^^i. 

The voltage-biased suspended nanowire structure 
shown in Fig. [9] was first analyzed in Ref. [41 for the 
case of nonresonant electronic transmission through the 
wire. Considering the electronic Fermi level to be far 
from the energy levels of the quantized longitudinal mo- 
tion of the electrons in the wire, the authors used pertur- 
bation theory with respect to the tunnel barrier trans- 
mission, Tg//, to calculate the current and conductance 



through the system. To second order in this expansion, 
the effective Hamiltonian coupling the electrons in the 
two leads through the virtual energy level of the wire can 
be shown to be of the form, 



a,k 



(16) 



k,k' 



Here, electrons in the leads are described by the first 
term while the fundamental mode of the fiexural vibra- 
tions^ is described by its harmonic oscillator Hamilto- 
nian (the second term in Eq. (16)). The additional phase 
factor e^^^^ connected with the tunneling electrons 
(third term in (16)) describes the magnetic field depen- 
dent phase which the electrons acquire between the two 
sequential tunneling events at the ends of the wire. The 
specifics of this term for the problem considered is that 
the electronic phase is not a c-number, but an operator 
acting in the quantum vibrational space. 

As the effective coupling strength for the tunneling 
electrons experiences quantum fiuctuations so does the 
number of different tunneling channels. Since these chan- 
nels obey quantum nanomechanical dynamics, we obtain 
a typical picture of entanglement between the quantum 
nanomechanical and electronic degrees of freedom. The 
outcome of this is that new polaronic states (referred to 
as "swinging states" in Ref. [411) are formed in the wire 
as a result of this entanglement, and charge transport 
can now be viewed as electronic transmission through 
these intermediate states. It is interesting to note that 
the strength of such polaronic coupling is determined by 
magnetic fiux and can thus be tuned by the external mag- 
netic field. 



From the effective tunneling Hamiltonian, Eq. (16), the 
electric current over the voltage-biased leads can easily 
be calculated. What is new in this analysis is that the 
matrix elements corresponding to electronic tunneling are 
now described through operators in the vibrational space 
of the wire, something that needs to be considered when 
tracing out the electronic and vibrational phase space. 
The resulting expression for the current is presented in 
Eq. m, 



oo oo 



n=Oi=-n 

X I de[/Ke)(l - fr{e - Ihw)) - /,(e)(l - fi{e - ^M)] • 

Here, Go is the zero field conductance, e is the elec- 
tric charge and the oscillating nanowire is assumed to be 



^ For simplicity, only the fundamental mode is considered in this 
analysis. 
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FIG. 10: Sketch of the different transmission channels avail- 
able for electrons tunneling through the oscillating nanowire 
of Fig. [9] Electrons with energy s tunneling from the left 
(source) to the right (drain) lead are transmitted in both elas- 
tic and inelastic tunneling channels (top image) with the cor- 
responding energy exchange, e + ihco; ^ = 0, ±1, ±2, — Due 
to Pauli-principle restrictions, some of the inelastic channels 
are affected by the electronic population in the drain lead 
(shaded region, lower image), which together with a reduc- 
tion of the tunneling rate in the elastic channel i = (see 
text and caption to Fig. [9]) leads to a finite magnetoresistance 
of the wire. The Pauli-principle restrictions are important 
only if e is close to the chemical potential in the drain lead. 
This is why the total current reduction saturates and becomes 
independent of both temperature and bias voltage for large 
enough bias voltage (see text). 



in thermal equilibrium with the environment, described 
through P(n) = (1 - e-^^^)e-^^^^, the probability of 
finding the oscillator in a quantum state n with energy 
nhu at a temperature T {/3 = {kBT)~^ and uo is the 
frequency of oscillation of the nanowire). The electron- 
vibron coupling is expressed through the coupling con- 
stant ^, which scales linearly with the magnetic field. 
Finally fi^r are the Fermi functions for the left and right 
lead kept at chemical potential jii^r = respectively 
with V the bias voltage. 

The above expression represents contributions to the 
current from a number of different inelastic tunneling 
channels, ^ 7^ 0, together with the elastic channel. Ana- 
lyzing Eq. (17) one finds that the destructive Aharonov- 



Bohm interference to the elastic channel caused by the 
quantum fluctuations in the magnetic flux (proportional 
to the area enclosed by the vibrating quantum wire) is 
compensated by the additional inelastic tunneling chan- 
nels. These inelastic channels increase the electronic 
transport through the system, and hence compensates 
the effect of the quantum suppression in the elastic chan- 
nel. The influence of these two effects on the total trans- 
port depends on the extent to which the inelastic chan- 
nels are permitted by the Pauli principle, which limits the 
electron-vibron energy transfer. The corresponding en- 



ergy diagram for the different electronic tunneling chan- 
nels is presented in Fig. [lO] 

Neglecting with the Fermi statistics of the electrons 
(high temperature limit), the restrictions imposed by the 
summation over the possible inelastic tunneling channels 
in Eq. ( 17) disappears, and the expression for the conduc- 



tance exactly coincides with the transmission through the 
non- vibrating wire (see Eq. ([l8|). In the low temperature 
limit however, the result of the above mentioned quan- 
tum mechanical suppression of the elastic channel and 
the Pauli restrictions on some of the inelastic channels 
is the appearance of a finite magnetoconductance. As 
a result, the strongest effects of the predicted quantum 
magnetoconductance occurs in the low voltage, low tem- 
perature limit: hco ^ eV^kBT. At higher energies, the 
Pauli restrictions become less important, and these effect 
only give corrections to the conductance as compared to 
the non- vibrating wire. Mathematically the asymptotic 
limits to the magnetoconductance are found as. 



_G 
Go 



(3hj > 1 

6^ (3hj<^i, 



(18) 



which can also be seen in Fig. ll(a)| Similar results 
of quantum mechanical effects on the conductance have 
been reported for the case of ballistic electron transport 
through a carbon nanotube encapsulating a single mov- 
able fullerene scatterei^^ (see also Ref. [43] for a resonant 
tunneling treatment of a similar system). 

The effects of the quantum suppression of the low en- 
ergy tunneling processes outlined above are only visible in 
the low temperature, low bias voltage limit and are hence 
hard to verify experimentally. If, however, the electrical 
current is instead considered these effects are indeed vis- 
ible also at higher temperatures and voltages. This can 
be understood as the Pauli restrictions on the inelastic 
tunneling channels affect only the low energy electrons 
(low bias voltage). Thus, the low voltage current is re- 
duced from the non- vibrating ohmic current, / = GqV ^ 
by an amount that is given by the extent to which the 
elastic channel is suppressed (which in turn depends on 
the magnetic field). As the bias voltage is increased the 
inelastic tunneling channels are opened and the current 
increases accordingly. However, a further increase in the 
bias voltage when the conditions eV ^ ^hx)^ ksT are ful- 
filled will not be affected by the Pauli restrictions due 
to the large energy scales of the electrons. As a result, 
these effects are displayed as a temperature- and bias 
voltage-independent current deficit, as compared to the 
non- vibrating wire, as recently reported in Ref. [40], 



lim / 



(19) 



The effect of the quadratic magnetic field dependence 
on the current deficit can be experimentally observed by 
extrapolating the high voltage asy mptotes of the I — V 
curves to the V ^ limit (see Fig. |ll(b)' ). For realistic 
parameters, we estimate that these effects should be 
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compared to the excitation strength set by the electronic 
emission of vibrons. If this is not the case, such electron- 
vibron coupHng cannot be ignored, even if the electronic 
tunneling is small. In order to account for the mutual 
electron- vibron dynamics, one should evaluate the evo- 
lution of the density matrix operating on both the vi- 
bronic and electronic subspaces, and, using this density 
matrix, calculate the vibronic contribution to the current 
as well as the electron- vibron contribution to the Lorentz 
force. This analysis was recently performed in Ref. [40] . 
where it was shown that the stationary density matrix 
can be found as a result of a balance between the elec- 
tron assisted emission and absorption of vibrons. The 
equation for the stationary reduced density matrix, p, of 
the electron- vibron subspace takes the form. 



[Hoscjp] = \Teff\ 



(ii + i2)p + p(j| + 4) 



(20) 
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■iXX 







012345678 

V [a.u.] 

(b) 

FIG. 11: (Color online) (a) Magnetoconductance, G/Go, 
through suspended SWNT system as a function of tempera- 
ture G = kBT/huj and magnetic field H/Hq. Reprinted with 
permission from [41], R. I. Shekhter et al, Phys. Rev. Lett, 
97, 156801 (2006). © 2006, American Physical Society (b) 
Predicted offset current as a function of bias voltage for two 
different temperatures and two different vibrational frequen- 
cies. Note that the offset current (red line) does not extrap- 
olate to the origin and scales as the square of the magnetic 
field, (j)^ oc . Data shown in arbitrary units. Reprinted 
with permission from [lO], G. Sonne. 



observable for bias voltages V ~ 30 pV and currents 
/ ^ 3pA in magnetic fields of ~ 20T. 



By exploiting the equilibrium distribution of the flex- 
ural vibrations, P{n) in Eq. (17), we have neglected with 
the non-equilibrium effects which can be stimulated in 
the vibronic subspace due to the coupling to the cur- 
rent carrying electrons. This assumption is valid if the 
internal relaxation of the vibronic subsystem is strong 



Here, Ji^2 are operators that take into account the 
electron-vibron coupling and Hq is the first two terms 
in Eq. (16). Also, the average Lorentz force (21a) and 



momentum (21b) on the wire can be found by multiply- 



ing Eq. ( 20 ) with the deflection and momentum operators 



and tracing out the nanowire's degrees of freedom, 

^{x) = \T,ff\'Tv ((Ji - J2)p + p{j\ - J|)) (21a) 

{p)=0. (21b) 

In Ref. [40] it was shown that although each inelastic 
electronic tunneling channel (see Eq. ( [Tt] )) is significantly 
renormalized by the non-equilibrium of the vibronic sub- 
system, the total high voltage limit, eV ^ kBT^hco^ to 
the current is still given by Eq. (19). Furthermore, any 



corrections to this expression from the thermal environ- 
ment was shown to decay exponentially in this limit for 
all relaxation strengths in the vibronic subsystem. 



B. Nanoelectromechanics of suspended nanowire 
superconducting weak link 

Superconducting ordering in the electronic subsystem 
qualitatively changes the electromechanical coupling in 
the suspended nanowire based NEM tunneling devices 
considered so far. If the electronic subsystem of the leads 
in Fig. [9] undergoes a superconducting phase transition 
a new type of electronic coupling between the leads oc- 
curs. This coupling allows for the possibility of coherent 
tunneling of a pair of electrons (a Cooper pair), forming 
a superconducting condensate of electrons in each of the 
two leads. As the process of Cooper pair tunneling is a 
dissipationless ground state property of the superconduc- 
tors, a finite dc current over the leads will be present even 
at zero bias voltage (the dc Josephson effect). If instead 
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a finite dc bias is applied an alternating Joseptison cur- 
rent is set up over the leads, representing the tunneling 
response of the junction to the bias, whose zero average 
value guarantees that no electric power is absorbed from 
the voltage source. The frequency of these current oscil- 
lations is set by the energy gained by a single Cooper pair 
during its transition between the superconducting leads: 
hu = 2eV (the ac Josephson effect). 

If now the tunneling Cooper pairs are coupled to the 
mechanical vibrations of the suspended nanowire, an 
electromotive force is induced due to the presence of the 
external magnetic field causing the vibrations^. The work 
done by this force on a tunneling Cooper pair renormal- 
izes the total energy gained by it and hence changes the 
frequency of the Josephson oscillations. On the other 
hand, the mechanical vibrations of the nanowire are af- 
fected by the Lorentz force which oscillates in time due 
to the ac variations of the Josephson current in the wire. 
From these considerations a set of equations determin- 
ing the Josephson frequency renormalization and the me- 
chanical motion of the wire were derived in Ref. [44] . 
These equations were found by considering the Andreev 
level formation in the NEM-superconducting weak link 
in the low tunnel barrier transparency limit (see also Ap- 
pendix [A| and are fully consistent with the qualitative 
picture describe above. Using the dimensionless coor- 
dinates y, e and 7 for the wire's deflection coordinate, 
driving force and phenomenological damping respectively 
these equation are. 



F + 7F + r = esm{(p) 
cp = V-Y. 



(22a) 
(22b) 



In (22), V is the dimensionless applied bias voltage and 
the driving force e ex H'^. Eq. 



([22| is thus the equation 
of motion for the deflection coordinate of the wire which 
is driven by a force proportional to a term reminiscent 
of the ac Josephson current, cx sin(Vt). Note however, 
that as the wire moves in the magnetic field it induces 
an electromotive force proportional to the time rate of 
change of the deflection coordinate, which counteracts 
the motion causing the vibration. This term is given 
by the second argument in Eq. ( |22b| ) (see Ref. [44 for 
details), and is responsible for the characteristic resonant 
phenomena described below. 



Numerical simulations of Eq. (22) shows that for 



small driving forces (low magnetic fields) the system will 
achieve resonant conditions only at the eigenfrequ ency of 
the fundamental mode of the wire F = 1, see Fig. 12(a) 



For larger driving forces however, also pa rametric exci- 
tations V = 2 can be achieved. Fig. 12(b) Furthermore, 



the amplitude of these modes is initially seen to be an 



^ The inductive coupling between flexural wire vibrations and the 
supercurrent flow in a SQUID loop was considered in Refs. [541 
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FIG. 12: Time- averaged vibrational amplitude of the super- 
conducting suspended nanotube as a function of driving volt- 
age, V, clearly showing the onset of the parametric resonance 
at higher driving force, e. Reprinted with permission from 
[44 , G. Sonne et al, Phys. Rev. B, 78, 144501 (2008). © 
2008, American Physical Society. 



increasing function of the magnetic field strength, but 
even tually saturates at some critical driving force e*. 
Fig. |13(a)[ To explain these phenomena we derive two 
stability equations for the amplitude and phase of the vi- 
brating nanowire which well describe the onset and sat- 
uration of finite vibrations at resonant conditions. These 
equations also show that for the system considered there 
exists a window of bistability in the vibrational ampli- 
tude around the resonance peaks. 

Being that the voltage source pumps energy into the 
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FIG. 13: (Color online) (a) Time- averaged vibrational ampli- 
tude as a function of driving force for the second resonance 
peak, V — 2^ clearly showing the saturation of the vibrational 
amplitude at e* ~ 0.005. Inset shows corresponding dimen- 
sionless dc current as a function of the magnetic field, (b) 
Predicted current bistability and hysteresis for two voltage 
pulses close to the parametric resonance. As shown, applying 
either a positive or negative voltage pulse in time can shift the 
system between the two stability points (1 and 2) which are 
separated by a finite current difference. Reprinted with per- 
mission from 44 , G. Sonne et al, Phys. Rev. B, 78, 144501 
(2008). © 2008, American Physical Society. 



the system considered should exhibit both negative and 
positive magnetoresistance as the current is first an in- 
creasing function of the magnetic field (increasing vibra- 
tional amplitude) but falls off as jdc ^ once the 
driving force reaches the critical value e* which indicates 
the onset of the bistability in the vibrational amplitude 
(see inset Fig. |13(a)] ). 

Also, we have shown that the bistability region around 
the resonance peaks can be utilized to directly probe the 
underlying quantum mechanics of the system. As was 
mentioned above, close to resonance conditions the sys- 
tem exhibits a mechanical bistability, i.e. two dynami- 
cally stable solutions of oscillation can be found for the 
same bias voltage but with different vibrational ampli- 
tudes. Which of the two stability point the system will 
correspond to will in general depend on initial conditions, 
however, by applying finite voltage pulses in time one can 
move the system between the two. Consider for example 
the situation shown in Fig. 1 13(b)] where the system is ini- 
tially found at the stable point 1, corresponding to the 
lower vibrational amplitude (lower current) and the bias 
voltage is slightly off the parametric resonance. Apply- 
ing two voltage pulses in time, (i) and (ii), will move the 
system along the phase space trajectory of the stability 
point defining it. As the region of bistability is only found 
within a small window of bias voltage around the reso- 
nance peaks, pulse (i) will in this example force the sys- 
tem into the second stability point, which will return the 
system to the higher stability point, 2, on the up-sweep. 
As the dc current scales directly with the vibrational am- 
plitude we thus predict that this pulse will result in a 
measurable current difference to pulse (ii) which will re- 
turn the system to the initial stability point. Considering 
realistic parameters we predict that the measurable cur- 
rent difference (1 ^ 2) ~ 5nA should be experimentally 
observable for magnetic fields H ^ 20 mT and bias volt- 
ages V ^ 5 /iV. 



IV. CONCLUSIONS 



oscillating nanomechanical system, it can be shown by 
direct integration of the equation of motion that at reso- 
nant conditions a finite dc current, jdc^ is set up over the 
system. This comes about as the mechanical energy of 
the nanowire is a constant of time at resonant conditions, 
hence any pumping of the system at these conditions will 
lead to an energy transfer through it, proportional to the 
real damping coefficient 7. As such, the resonant dc cur- 
rent can be shown to scale as, 



Jdc CX 



(23) 



where time-average rate of change of the deflection coor- 
dinate, (y(t)^), behaves phenomenologically equivalent 
to the deflection coordinate. As such, we predict that 



Coupling between electrical and mechanical degrees of 
freedom is the basic mechanism behind the functional- 
ity of any nanoelectromechanical system. Such coupling 
can be achieved in many different ways. In this Re- 
view we have focused on the effects of varying the spatial 
concentration of the electronic charge or current in sys- 
tems incorporating suspended nanowires as the mechan- 
ical element. These structures are particularly interest- 
ing to study as they can be used as efficient electrome- 
chanical transducers that serve simultaneously as electric 
weak links for the tunneling electrons and as nanome- 
chanical resonators. This allows one to use the mutual 
coupling between the electric current oscillations in time 
and the mechanical vibrations to achieve a highly non- 
linear electromechanical coupling if strong driving volt- 
ages or magnetic fields are applied. The resulting self- 
supported nanoelectromechanical dynamics can be main- 
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tained through the development of an electromechanical 
instability, caused either by the shuttle phenomenon (for 
systems with varying charge concentration) or through 
resonant pumping of nanovibrations by the electric cur- 
rent flow (in systems with strong current concentration). 

In many of the systems considered in this Review, 
quantum coherence is an essential feature of the discussed 
nanoelectromechanical operations. For example, we have 
shown that quantum effects in the electronic subsystem, 
such as coherent electron transmission through a double 
barrier structure, can be coupled to the quantum coher- 
ence of the mechanical flexural vibrations (Section |lII A ). 
The resulting electromechanical entanglement establishes 
"swinging" polaronic states for the tunneling electrons, 
which allows the quantum vibrations to be investigated 
by interferometry. Also, we have shown that the exis- 
tence of more than one mode of nanomechanical bending 
of the suspended nanowire, a direct consequence of the 
spatial extension of the mechanically vibrating element, 
makes the operation of the NEM device even richer. As 
an example we mention the multimode shuttle instabil- 



ity discussed in Section II D| where it was shown that by 
applying a "local" external probe in the form of current 
injection from an STM tip the "global" mechanical prop- 
erties of the system as a whole can be altered. In par- 
ticular, the competition between the different mechanical 
modes of vibrations will lead to a self-organization of the 
shuttle vibrations of the suspended nanowire. 

Experimental realization of the phenomena dis- 
cussed in this Review is an experimentally challenging 
task. Although NEM-structures based on suspended 
nanowir es have been m anufactured by a number of 

l^bsl3 14 16|17 | 18|19 | 22|4 5l46r47 48^ ^^^^ fo^ g^^^^g 

nal magnetic fields and applied voltages offers a real ex- 
perimental challenge. Another obstacle to experimental 
realization of these systems is damping of the nanome- 
chanical vibrations. If, for example, the mechanical 
damping is too large, then the voltage threshold required 
for the onset of the electromechanical instability may well 
exceed the thermal limit, set by Joule heating of the sys- 
tem, and the system might burrP^. However, for the 
NEM systems considered here typical experimental pa- 
rameters such as quality factors of the order Q ~ 1000 
and bias voltages of order 1 V the suggested multimode 
shuttling and superconducting pumping of nanovibra- 
tions should be experimentally observable. The effects 
of the interferometry of the quantum nanovibrations dis- 
cussed here do nevertheless demand the presently best 
achievable experimental conditions. 

In this Review we have discussed the coupling of meso- 
scopic electrons to nanomechanical radio frequency vi- 
brations. It is interesting to compare this to how meso- 
scopic electrons couple to external electromagnetic fields 
of much higher frequencies. Such coupling has been 
been shown to result in a number of interesting non- 
equilibrium mesosc opic p henomena in both normat^^^ 
and superconducting^^ mesoscopic structures if the mi- 
crowave photon energy is comparable to the electronic 



energy-level spacing. The much lower mechanical fre- 
quencies of NEMS devices, on the other hand, can easily 
be tuned to match the smaller electronic energy scales 
that characterize the widths (due to tunneling) of elec- 
tronic energy levels, the applied bias voltage, or the 
frequency of Rabi oscillations in the population of en- 
ergy levels. The fact that electrons in mesoscopic sys- 
tems can couple to oscillations corresponding to such 
very different energy scales, makes it tempting to spec- 
ulate that mesoscopic electrons might be used for mak- 
ing microwave-mechanical transducers by simultaneously 
coupling them to both electromagnetic and nanomechan- 
ical degrees of freedom. If such an electron-mediated 
microwave- mechanical coupling is possible in principle, 
one may wish to consider quantum coherence effects in 
structures where both electromagnetic, mechanical and 
electronic degrees of freedom have been quantum me- 
chanically entangled. The extent to which such specu- 
lations are realistic is an open question which we leave 
for future research in the exciting field of nanoelectrome- 
chanics to answer. 
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APPENDIX A: DERIVATION OF THE 
EXPRESSION FOR THE FORCE ON THE 
PUMPED NANOMECHANICAL SYSTEM 

Consider the Hamiltonian, for the electronic sub- 
system of the pumped Josephson vibrations of the sus- 
pended nanotube, Eq. (1) in Ref. [44] . 

n = J dx¥{x) (Ho + Ha) ^(^) 



Ho 



2m 



d_ 
dx 



ieHu{x^ t) 

' h 



-CT,U{X) 



Ha = A(x) [ax cos (j){t) + sgn(x)cr^ sin ^(t)] . (Al) 



In ( Al ), ^"f" [^] are two-component Nambu spinors, cr^ are 
the Pauli matrices in Nambu space and the deflection of 
the tube is given by u{x^t) = u{x)a{t)^ where u{x) is 
the dimensionless, normalized profile of the fundamen- 
tal bending mode and a{t) determines the amplitude of 
vibration. The potential U{x) describes the barrier be- 
tween the nanotube and the bulk superconducting elec- 
trodes where the gap parameter is A(x) = Ao6(2|x| — L) 
and (j){t) — 2eVt/h is the phase difference across the 
junction. As indicated in Ref. [44 a gauge transform, 
^iSfi^-iS ^-^^Yi S = cr^eH u{x' ,t)dx' /h, shifts the de- 
pendence of the vector potential induced by the nan- 
otube deflections from the kinetic part of the Hamiltonian 
to the phase differences over the leads, (j){t) (p{t) = 

(j){t) - a(t)4.eH J^^^ u{x)dx/h. 
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Suppose now that we have an eigenstate of the system, 
IV^(t)), such that IV^(t)) = e-^/o^o(^')d^'/^(|^o(t)) + |^(t))) 
where \ipo{t)) is the ground state of the system corre- 
sponding to energy Eo{t) and is some state or- 
thogonal to the ground state. To be able to evaluate 
the force on the system with respect to the ground state 
we need to show that is small, which can be done 
by expanding it in a complete set of eigenstates of the 
Hamiltonian, \i^{t)) = An{t)\ipn{t)) . Inserting this 
form into the Schrodinger equation and multiplying from 
the left with the state {ipnit)), we find that the ampli- 
tudes An{t) can be expressed as. 



An{t) = 



En{t)-Eo{t) 



(A2) 



To evaluate this further we need to find and expression 
for the term in the nominator of Eq. ( A2 ) which can be 



done by considering the time derivative of the eigenstate 
solution Hltpoit)) = Eo\ipo{t)) and evaluating the matrix 
elements of this equation with Performing this 

analysis we find that we can express Eq. (A2) as. 



An{t) 



ih{Mt)\{dtn)Mt)) 

{Er,{t) - Eo{t))^ 



(A3) 



As the only time-dependence of the transformed Hamil- 
tonian is found in the phase difference over the leads we 
identify that the nominator of this expression is propor- 
tional to the bias voltage, dtH cx eV^ whereas the denom- 
inator is the spacing of the energy levels in the systems, 
(X Aq. Hence, for the system considered eV <C Aq, we 
can safely evaluate the expression for the force with re- 
spect to the fixed-phase ground state |?/^o(^))- 



The force on the wire is then found by evaluating 
the expectation value of the force operator F cx dH/da 
which can easily be shown to be F = (^/^ol^lV^o) = 
—dEo[(p{a)]/da, where Eo{(p) is the ground state energy 
as given in Ref. [44 . Similarly the Josephson current can 
be shown to be j = {2e/h)dEo{(p)/d(p and we recover 
the total expression for the force on the oscillating wire 
as given in Eq. (22). 
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